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ABSTRACT 

A microwave radar technique for remotely measuring the vector wave 
number spectrum of the ocean surface is described. The technique, 
which employs short -pulse, noncoherent radars in a conical scan mode 
near vertical incidence, is shown to be suitable for both aircraft 
and satellite application. The technique has been validated at 10 km 
aircraft altitude, where we have found excellent agreement between 
buoy and radar-inferred absolute wave height spectra. 
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1. INTRODUCTION 

For several years now, we have been endeavoring to develop a 
microwave radar technique for measuring ocean wave directional spectra 
that would be suitable for satellite application. Basically, we have 
been seeking to define an alternative to the coherent imaging radar 
approach that wa2 adopted for Seasat, the nation’s first oceanographic 
satellite (Be^ et aK, 1981). Our motivation has been to find an 
alternative measurement approach that would at the same time a) be 
simpler and less costly, b) be capable of truly global measurements, 
and c) be more accurate. 

In this we believe we have been successful . Theoretically, and on the 
basis of aircraft flight experiments we have determined that such 
global -scale satellite measurements are feasible. The measurements 
can be made with relatively simple, noncoherent short-pulse radars 
operating In a conical scan mode near vertical incidence, Q ~ 10°. 

No new technological developments are required. Rather, these 
measurements can be made with existing space-qualified hardware. 

For example, with some relatively minor modifications such as the 
addition of a modest-gain scanning antenna, the Seasat altimeter can 
be adapted to perform these measurements. The measurements are 
inherently of high resolution spectrally in both wave number and 
direction, and as we shall see, they will be remarkably accurate as 
well . 

A typical satellite measurement geometry is illustrated in figure 1. 
For the assumed satellite altitude of 700 km and Incidence angle of 
10°, the radius of the scan pattern on the ocean surface is 
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approximately 130 km. A 3-rpm antenna rotation rate is selected as 
a reasonable compromise between coverage and Integration time 
requirements. The measurement cells (not to be confused with the 
instantaneous field of view, or antenna 'footprint') are roughly 130 
km squares situated one on either side of the subsatellite track. 
Basically, the measurement product consists of two statistically 
stable estimates of the polar-symmetric vector wave number spectrum, 
one on either side of the subsatellite track. If less than 180° of 
look is allowed, then these measurements can be confined to an area 
considerably smaller than the nominal 130 km square as is evident 
from figure 1. 

Although the technique we shall be considering employs short-pulse 
waveforms, it is not in its most fundamental aspect different from 
the two-frequency technique investigated theoretically by Alpers 
and Hasselmann (1978) and experimentally by Johnson et a]_. (1981). 

In both techniques the basic measurement principle is the same. 

This is the directional selectivity that results as a natural 
consequence of the phase-front matching of electromagnetic and ocean 
wave components. The choice of waveforms, and the manner of detection, 
is however a critical one. Jackson (1981)— hereinafter referred to 
as J— has shown that the narrowband two-frequency technique has. 
Inherently, a very low signal -to-noise ratio (SNR) compared to the 
short-pulse technique. Basically, this is because the sea-spectrum 
is relatively broadband whereas the two-frequency beat -wave signal 
is comparatively narrowband. For large footprint dimensions this 
results in modulation signal power being detected only in a very 
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narrow spectral band, and consequently, the signal energy Is small 
compared to the fading variance. 

Our work differs from that of Alp ere and Hasselmann (1978) In another 
Important respect. This Is In the choice of incidence angles. 

Alpers and Hasselmann (1978) were concerned with large-angle 
measurements whereas our concern is with small angles of Incidence. 
There are several reasons why we have chosen to study small -angle 
scatter. First, as should be apparent from the above discussion of 
the measurement geometry, small angles of Incidence are necessary 
at satellite altitudes in order to keep the scan radius to a minimum. 
If the nadir angle is too large, the scan pattern on the surface may 
exceed the scale of homogeneity of the wave field. Second, the 
reflectivity modulation mechanism in near-vertical backscatter is 
simpler and more predictable than it is in large-angle backscatter. 

In the near-vertical, specular backscatter regime the contrast 
modulation does not depend on the strong— and essentially 
unpredictable— hydrodynamic modulation of the short Bragg-diffracting 
water wave. The modulation mechanism is primarily a geometrical 
tilting effect, and consequently, it is more amenable to accurate 
modelling. Another reason for choosing small incidence angles is an 
obvious one, that the greater cross-section and lower link loss near 
nadir demands less transmitter power and antenna gain. This is an 
important consideration in the wide-band measurement approach that 
we are advocating. 

In the following we will discuss the three major conceptual elements 
that constitute the measurement technique, namely, 1) the principle 
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of directional selectivity, 2) the modulation mechanism in near- 
vertical backscatter, and 3) the use of short -pulse waveforms to 
detect the range reflectivity modulation. The discussion is Intended 
to provide a basic understanding of the measurement technique and to 
provide such results and formulas as will be found useful In the 
analysis of the aircraft data. For a fuller and more detailed 
theoretical treatment, the reader is referred to J, 

2. THE MEASUREMENT TECHNIQUE 

2.1 The principle of directional selectivity. 

We are concerned with fairly narrow antenna beams in a high-altitude 
measurement geometry. The relevant geometry is illustrated in figure 1. 
The situation desired is one where a) the antenna footprint is large 
compared to the scale of the waves, and b) the curvature of the 
wavefront is small compared to the directional spread of the waves. 

Now obviously, if the lateral beam spot dimension is large compared 
to the scale of the waves, then the waves cannot be resolved in 
azimuth (short of resorting to synthetic aperture). Rather, the wave 
contrasts will be averaged laterally across the beam. What is the 
effect of this lateral averaging? To understand the effect, imagine a 
Fourier decomposition of the two-dimensional reflectivity field 

into an angular spectrum of plane contrast waves. (The reflectivity 

# 

field can be imagined to be that measured by a very high resolution 
real -aperture imaging radar looking in the same azimuth direction.) 
Referring to figure 2, it is apparent that the effect of the lateral 
averaging is to eliminate or 'cancel out* any plane surface contrast 
wave that is not aligned with the beam direction. Only those surface 
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waves whose phase fronts are 'matched' to the electromagnetic (em) 
phase front can survive the lateral averaging. The effect of the 
broad footprint is then to isolate, or resolve surface contrast wave 
components whose wave vectors = (K, $) are aligned with the beam 
di recti on . 

The directional resolution is determined by and limited by a) 
the finiteness of th»? beam spot size in azimuth Ly, and b) the 
curvature of the wave front within the beam spot. If we assume; 
a Gaussian-shaped azimuth gain pattern, 

G(y) = exp(-y 2 /2Ly 2 ) (2.1) 

then it follows (e.g., from the Fresnel zone solution in 0) that 
the directional resolution 6$> , defined as the half-power spectral 
window width In azimuth, is given by 

6$ ~ 6Ky/K = a^TTnF [ ( KLy ) ” 2 + (L y cote/2H)2]i/2 (2.2) 

where H is the altitude. The first and second terms in (2.2) derive 
respectively from the finite-footprint and wave-front curvature effects. 
In our aricraft experiment geometry, H ~ 10 km, 6 ~ 13°, and Ly ~ 300 m. 
(half power width Ly* - 2 / 21 n2 Ly ~ 700 m) . For a typical 200 m water 
wave, we have <$$ ~ 17°. In a typical satellite measurement, H = 700 km, 

9 = 10°, and Ly = 8.5 km (L y * = 20 km), in which case 5$ ~ 5°.* 

2.2 The reflectivity modulation in near-vertical backscatter. 

Near vertical incidence, 0 <15°, microwave backscatter from the sea 
occurs by means of quasi -specular reflections from wave facets 


*The directional resolution quoted in J (equation 79) is wrong. 


6 


oriented normal to the radar's line of sight. The average backscatter 
cross -section o° Is proportional to the probability density function 
(pdf) of orthogonal surface slopes satisfying the specular condition 
for backscatter: 3c/3x * tane ; ss/9y a 0. The cross-section is given 
by (e.g., Valenzuela, 1978); 

a°(0,$) « pTrsec 4 0p(tane,O) (2.3) 

where p is the slope pdf expressed in the radar's coordinate system, 
x in the plane of Incidence, arid where p Is a diffraction-modified 
normal Incidence Fresnel reflectivity (Brown, 1978). 

H/drodynamic modulation is a second-order effect in near-vertical 
backscatter. Consider that, firstly, for most microwave frequencies, 
the most strongly forced waves, the gravity-capillary waves, lie 
under the diffraction limit (about three em wavelengths in the 
horizontal according to Brown (1978)). Thus, they are only weakly 
sensed, and to the extent that they are, it is via a diffuse 
diffraction field that can be only very weakly modulated by 
geometrical tilting. Secondly, the specular component derives from 
the entire wave ensemble, including waves on all scales, from the 
scale of the dominant waves we are seeking to measure down to the 
scale of the diffraction limit. For this large ensemble of waves, 
it is reasonable to assume that hydrodynamic forcing and wave-wave 
interaction effects are of secondary importance. To the extent that 
hydrodynamic nonlinearities effect the em modulation, they are to be 
attributed to the entire wave ensemble rather than a particular 
water wave component. Neglecting second-order effects, the surface 
can be treated as a free-wave superposition possessing Gaussian 
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statistics. If the large-wave slopes are then assumed to be small 
compared to the total rms surface slope, the modulation can be 
modelled by the following linear 4 t11t model'. 

The backscatter cross-section of a small patch of sea surface of 
area A (cf. figure 3) Is given by o « o°A, where the normalized 
cross-section a 0 is assumed to be the average a 0 of the sea surface 
in a tilted reference frame. Thus, if 0' and are the local 
Incidence and azimuth angles, we suppose that a°(patch) - o°(e',0'). 
For small large-wave tilts 6, the fractional cross-section variation 
is given by 


6a a 6a° + 6A (2.4) 

"a ~o° 

The elementary surface area is that area contained in the range 
Interval cAt/2, To first order in 5, A is given by A « Ay(cAx/2)csc 0' . 
Provided that 6 « 0, the local Incidence angle can be approximated 
by 0' ~ e - 8?/3x. Thus, to first order in 5 it follows that 
6A/A = cot0 8c/8x. Since the azimuthal dependence of a 0 is 
small compared to the 0 dependence, it follows that the tilt term 
6a°/a° is also proportional to the large-wave slope component 
in the plane of incidence. From (2.3), 

6a° a __ JL 3fi 8c + 0(6 2 ) (2.5) 

a 0 P stane 3x 

The fractional range -reflectivity modulation seen by the’ radar is 
6t >; o averaged laterally across the beam: 
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m(x,$) * 


y (y )dy 


( 2 . 6 ) 


The directional modulation spectrum Is defined by 


P m (K,«) - (2rr)“ l /<m(x,$)m(x + K ,*)>exp(-1Ke)<fe (2.7) 

where the angle brackets denote ensemble average. Now le* G b& given 
by the Gaussian pattern (2.1), and consider the limiting case of 
very large footprints, KLy » 1. It is easy to show then that P m 
is proportional to the directional slope spectrum as 


MM) 


rs 

~i — 

u y 


fcote " ilHfiJ 2 K Z F(M) 
3 tans J 


( 2 . 8 ) 


where F is the two-sided, polar-symmetric height spectrum defined 
so that the height variance, 


<C 2 > s /o7o*2F<M)KdKd* (2.9) 

The rms modulation depth, by definition, is given by 

u(») * <m 2 (x,.)> w2 » [/„”2P m (K,«.)dK]i' 2 (2.10) 

t 

It should be pointed out that, strictly, the large footprint 
limiting form is valid only if Ly is much larger than the lateral 
decorrelation scale of 3c/3x. This is equivalent to the condition 
KLy >> 1 in the general case of directionally spread seas, but not 

in the case of unidirectional, long-crested swell. We have encount- 

** 

ered such a swell in our aircraft experiment, where the crest length 
was very long compared to the antenna beamwidth. In such a case a 
separate calculation must be carried out, one which accounts for the 
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curvature of tne em wave front* But as the case we encountered was 
exceptional among our data, we have neglected to perform such a 
calculation. 

Now It Is only consistent at this point to assume that the slope 
pdf Is Gaussian. Indeed, It would be Inconsistent to assume other- 
wise since the tilt model Is predicated on an assumption of free, 
non-interacting waves, and this can only Imply normal statistics. 

The K u -band scatterometer data of Jones et aj_. (1977) analyzed by 
Wentz (1977) show in fact that the pdf Is nearly normal. More Inter- 
esting though, the data Indicate that the pdf at K u -band frequencies 
is very nearly Isotropic. This Is convenient, as It simplifies the 
measurement of bl- or multi-modal directional spectra since the 
sensitivity Is independent of azimuth, and no relative weighting of 
different directional components is required. If the slope pdf Is 
Gaussian and Isotropic, then the sensitivity coefficient, the factor 
of K?F in (2.8), can be written as 

a 5 /TFrcote +2tan0T 2 (2.11) 

17 l <JWF» 

where <|vc| 2 > Is the mean square wave slope effective at the 
particular radar operatng frequency (diffraction-effective mean 
square slope). 

The linear tilt model solution (2.8) is Identical to the first term 

i 

in the series expansion of the geometrical optics solution obtained 
by J. The second-order terms consist of a em and a hydrodynamic 
(hydro) term. The two terms are of comparable magnitude, both'scaling 
as the large-wave steepness < 5 0 = Kq ^ 2 ) 172 to the fourth power. The 
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em term Is Independent of hydrodynamic nonlinearity and arises In 

i 

scattering from a normally distributed sea surface. The hydro term 
Is due to the non-Gausslan statistics associated with hydrodynamic 
nonlinearity, and Is given In terms of various third moment statistics 
In wave height and slope. Since these statistics scale with S 0 
(see J; also Huang and Long, 1980), the result Is that both em and 
hydro terms scale as 6 0 4 . The calculations of the second-order 
em term carried out in J Indicate that, firstly, the term Is 
generally small, and secondly, that the least harmonic distortion 
occurs In the neighborhood of 10° Incidence. The smallness of the 


second-order 


terms requires 


i.L.1 A. I J ,..l JJ.J , _ I 1 J 

mat one luiiuwiny umCjUai ivmS snuUiu 


be satisfied: 


1) <$ 0 cote « 1 

11) 6 0 tane « 1 

<\vW> 


( 2 . 12 ) 


If 1) Is violated seriously, an obvious consequence Is that the 
phase-front, or pulse, may Intersect the surface at more than one 
point. A less extreme but more general consequence of violating 
1) Is the confounding of the surface range coordinate with the wave 
height. The range coordinate will suffer a displacement Sx = ?cote. 
The net result will be a dispersion of the range coordinate by an 
amount <£ 2 > 1/z cote. This dispersion will represent a limit to 
the smallest wavelengths observable by this technique. Since 
< ? 2 >i /2 - 6 0 /« 0 , it follows that the upper limit on wave number as 
a function of the peak wave number is of the order of 
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w<0 “ («o cot0 )‘ 1 (2.13) 

For example, if 0 « 10° andfi 0 = 0.05 (fully aroused seas), then 
Kfnax ~ 3.5 Kq. In steep developing seas, ^ ~ 0.1, in which case 
Kmax ^ 1*75 1^. To illustrate the nature of the spurious response 
associated with the violation of 11), take the extreme case of swell 
under calm conditions. Obviously, if tane > S 0 then no backscatter 
occurs since there are no wave slopes satisfying the specular con^ 
dition. If tane < S 0 , the backscatter will now occur in periodically 
spaced bursts at points on the swell profile satisfying the specular 
condition. The backscatter will look like a string of delta functions, 
and will bear little resemblance to the swell profile save in its 
periodicity. Clearly, for the measurement to have decent fidelity 
(to the slope spectrum) there must be sufficient small-scale rough- 
ness, or in other words, a sufficient density of specular points. 
Practically, this means that the local wind speed should be in 
excess of several meters per second. 

Some guide to the selection of the 'best' incidence anqle may be had 
by the following. Assume that the inequalities (2.12) carry equal 
weight, that is, assume that the consequences of violating 1) and ii) 
are equally undesirable. Then we can minimize the sum i) + ii) with 
respect to 9 . This yields tane = <jvsj 2 >. For example, if the wind 
speed is 10 ms _i then using (4.8) we get e = 10°. More work alonq 
the lines established in J is required to get a better idea of what 
is really the best angle for minimizing the measurement nonlinearities. 
Unfortunately, the aircraft data are of little or no use to us here. 
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This Is because at the relatively low aircraft altitudes, the 
elevation beamwidth must be fairly broad in ordef’ to generate a 
sufficiently large beam spot for wave number resolution. In our 
aircraft experiment geometry, the 10° elevation beamwidth makes it 
virtually Impossible to establish the optimal angle since the likely 
range of 0 lies within the beamwidth. 

2.3 The short -pulse technique. 

In principle, the range reflectivity modulation spectrum P m (K ,<3> ) 
can be measured by either short-pulse or two-frequency techniques. 
However, as shown by J, the narrowband two-frequency technique has, 
inherently, a poor measurement SNR (signal -to-noise ratio) compared 
to the short -pulse technique. This is due to the use of narrowband 
waveforms that completely fill the beam. The analysis bandwidth <SK 
in this case is equal to the reciprocal of the range footprint 
dimension— the 'record length 1 . Hence the SNR « P m (K)6K will 
neccessarily be small when the footprint dimension is large. Since 
P m * Ly" 1 and SK « L^" 1 it follows that the two-frequency SNR is 
inversely proportional to the footprint area as noted by Alpers and 
Hasselmann (1978). 

In the short-pulse technique,* wide bandwidth, short pulses are 
used to resolve the wave structure in range. Backscattered pulses 
are integrated in surface-fixed range bins, and the range modulation 
spectrum is computed digitally from the observed sample of the 
range modulation m(x,$). For narrow pencil beams, the curvature 
of the wavefront can be neglected, and the surface range can be 
taken to be a linear function of the signal delay time t. If the 
*The terminology may be due to Tomiyasu (1971). 
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motion of the platform Is accounted for, then 

x = ct/(2sin0) + Vtcos ft s x 1 + Vtcost (2.14) 

where c is the speed of light, V is the platform speed, and where x 
and t are referred to the center of the beam spot at t “0, 5he 
coherency of the radiation results in random signal fading akin to 
the/ speckle observed when a coherent laser illuminates a 'rough' 
surface such as an ordinary piece of bond paper. The backscattered 
field statistics are complex Gaussian, the amplitude is Rayleigh- 
distributed, and the detected power is exponentially distributed 
(Moore et aU, 1975). If the measurement integration time is short 
{< 1 s), the surface can be regarded essentially as frozen. Consider 
the backscatter of a short pulse of length At. The surface range 
resolution cell is given by 

Ax e cAT/(2sin0) (2.15) 

It is assummed that Ax is small compared to the dominant wavelength. 

The backscattered power in a pulse transmitted at a time t^ can be 
modelled as a weakly modulated noise process of the form 

Wj(x 1 ,*) s W 0 (x',«t)[l + m(x,<t)]w(x,t 1 ) (2.16) 

where W Q * G 2 (x)a°(e,<S>) is the average backscattered power profile 
profile and w is the unit exponential fading process. The 
decorrelation length of the w-process is equal to the range 
resolution; more generally, if excess bandwidth is employed, the 
decorrelation length is given by the reciprocal of the pulse bandwidth 
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(Moore et_aj_., 1975). Since the autocorrelation function of the 
fluctuating component of the w-process Is concentrated near the 
origin, It follows that the spectrum of the (normalized) 
backscattered power Is given by 

Pl(M) - «(K) + P ln (M) + [1 + m 2 (4)]P w (K) (2.17) 

where <5 is the Dirac delta function. The fading spectrum is 
essentially the spectrum of the transmitted pulse. This follows 
since the backscattered field is given by the convolution of the 
pulse waveform with the surface impulse response, which, in the 
absence of modulation, is a complex Gaussian white noise process. 

For a Gaussian pulse shape with half power width At one finds (cf. 

J, eqs. 5 and 70 et seq .) : 

Pw(K) * Ax e xp[-(KAx) 2 /81n2] (2.18) 

2/2irTnT 

An integration of N independent pulses will reduce the fading variance 
by a factor of N“ l . Because of the platform motion, the backscattered 
pulses must be integrated in range bins that are fixed in the surface, 
otherwise the wave contrasts will be smeared out by range walk. This 
can be accomplished simply by delaying or advancing the trigger 
signal to the sampling gates according to the 1 1 ne-of-si ght relative 
speed between the platform and the surface.* The N pulse average can 
be expressed as 

N 

W n (t) = N -1][ Wj (t + ttj) (2.19) 

i=l 

*For satellites, it may prove difficult to specify the relative speed 
as a function of azimuth to the required accuracy. In this case, an 
alternative processing scheme, similar to that employed in the dual 
frequency technique, may have to be devised. 
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Where the rate of change of signal delay r a -(2V/c)sinecos$ . 
Omitting the dc term a.nd neglecting then the spectrum of the 
N-pulse average is given by 

P N (M) * P m (K,4>) + r%(K) 

If thermal noise is negligible, then the SNR is just the ratio of 
the signal spectrum to the residual fading spectrum, SNR “ NP m /P w . 
Using (2.18) for P w and assuming K « 2 tt/ax, we have 

SNR = 2/2717)? N P m (JK) 

AX 

The number of Independent pulses depends on the pulse repetition 
frequency (PRF), the Doppler bandwidth B^, and the integration 
time If the PRF > 2Bj, the signal is essentially continuously 

sampled and hence N = B^Tint. If the PRF « Bj, then the individual 
pulses are independent, in which case N = PRF X Tj n t. The Doppler 
bandwidth is determined by the interference rate of waves back- 
scattered from the lateral extremities of the range resolution cell 
From elementary considerations, or from equation 72 in J, 

Bd = (2V/X)e t#J ls1nO| 

Here B^ is the half power, post -detection Doppler spread in Hertz, 

X is the em wavelength, and s (Ly*/H)cos6 is the half-power 
azimuth beamwidth. 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 


The measurement integration time is limited by the azimuth scan rate 
which in turn is driven by coverage requi rements . The Integration 
time should not be longer than the time it takes to move one footprint 
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dimension. The modulation signal can only be built up coherently when 
the radar is viewing the same portion of the surface. When the beam 
moves to view a new, statistically independent patch of sea, the 
range modulation signal will evolve randomly, and further integration 
will proceed in an incoherent fashion, not only with respect to the 
scintillation or fading noise, but with respect to the modulation 
signal as well . Thus both P m and P w will be driven down as N" 1 . 

Thus, as the beam moves to view a new piece of the surface, the 
signal strength goes down as 1/Tj n t while the SNR approaches an 
asymptotic value. Since the antenna rotation is generally more rapid 
than the beam's translation, the azimuth scan rate determines the 
choice of integration time. Let us arbitrarily require that the beam 
move no more than one-half of its azimuth dimension. Then the 
integration time is set by 

Ti nt < A#/2* = fycsce/2* (2.23) 

An interesting consequence of (2.23) is that the SNR is independent 
of the footprint dimensions and hence of the antenna gain. This 
follows since the integration gain N * T-f n -t 05 Ly, while the signal 
spectrum P m « Ly" 1 . Thus, while the azimuth beamwidth affects 
the modulation signal strength (weakly as Ly" 1/2 ) it does not 
affect the measurement SNR. 

The number of degrees of freedom (DOF) in a measurement of P m (J<) 
is determined by the number of elementary wave number bands 
6K ~ 2ir/L x * contained in the spectral estimate. For example, 
consider an analysis with 25% resolution. Then the DOF of the 
estimate is given by (Blackmann and Tukey, 1953): 
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DOF - 2(0 ,25)/<5K ~ KL x */4tt (2.24) 

For example, if L x * = 20 km and K = 2tt/200 m then the DOF ~ 50, j 

3 

i 

i 

f 

3. A SATELLITE SYSTEM \ 

1 

These measurements can be made with a modified Seasat class ; 

radar altimeter. The pertinent Seasat altimeter characteristics 

are (Townsend, 1980): I 

Frequency: 13.5 GHz 

Pulse type: Linear FM, 1000:1 pulse 

compression 

Pulse length: 3.2 ns compressed 

Peak power: 2.0 KW 

PRF: 1000 Hz j 

Detection: Noncoherent square law j 

$ 

One can modify the Seasat altimeter* in such a way that it can 

perform a dual function, first as an altimeter per se , and second 

as a 'directional wave spectrometer 1 . In the conventional altimeter 

mode, mean altitude and wave height are determined from the delay 

time and broadening of the leading edge of the averaged return of 

nadir-directed pulses. There are several ways whereby transmitted 

pulses may be shared between the instrument's nadir altimeter mode j 

and off-nadir spectrometer mode;--for example, by power dividing or 

time sharing. Modification would entail the addition of a separate 

receiving section (post IF) and microprocessor as well as a separate 

rotating antenna. Pulse tracking, integration and spectral analysis 

Unfortunately, the Seasat engineering unit is no longer available j. 

for this purpose. j 
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functions would be Incorporated in the separate microprocessor. As an 
example, let us consider adding aim diameter, 3 rpm rotating antenna 
to the existing instrument. If we assume a 700 km satellite altitude 
and 10° nadir angle, then the measurement geometry Is that of figure 1, 
and the relevant measurement parameters are: 


Velocity : 

Beamwidth : 

Spot size; 

Rotation rate: 

Danno mertl 11 + An • 

nuiiHVi I wuy iuui wri • 

Doppler bandwidth : 
Integration time: 


V a 7 kms" 1 
Pe * B 1.6° 

L x * - L y * = 20 km 
(L x ~ Ly a 8.5 km) 

9 = 360720 s 

Av rs 9 .0 m ffmin 9_1K\ 

Bd « 18 | si n<3>| KHz (from 2.22) 
= 0.26 s (from 2.23) 


The PRF equals Bj at $ ~ 3° of foreward or aft. For most azimuths 
the PRF « Bd so that the number of independent samples is given by 


N = PRF X Ti n t = 260 = + 24 dB 

For illustrative purposes let us assume a Phillips* cutoff spectrum 
with a cos' 4 spreading factor: 

F(M) = U.005(4/3Tr)cos l4 («I»-<t 0 )K“ 1+ , K > K n 

l 0 , K < Kq 

Assume a 200 m water wavelength and upwave/downwave looks. Let 
the mean square slope as a function of windspeed be given by 
(4.8) and let U = 10 ms" 1 . Then we have (cf. 2.8, 2.11, and 2.21): 

P m = (2.95 X 10-* 4 m" l )(5.67 + 9.53) 2 (2.15 m 2 ) = 0.15 m 
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\i * 10 % 

SNR = 260 X 0.22 * +18 dB 

The directional resolution given by (2.2) is 6$ = 4.7°; but this 
assumes no rotation of the beam. Since the beam moves about 5° during 
the integration time, the actual resolution will be somewhat less. 

We have already calculated the DOF in the last section; DOF ~ 50. At 
10° incidence, o° ~ + 5 dB and is very nearly independent of windspeed. 
A link equation assuming 3 dB in losses and a noise factor of F = 6 
dB gives a signal -to-thermal -noise ratio of + 6 dB. Thus thermal 
noise is not a problem, even if half the transmitter power is shared 
with the altimeter mode. 

The spectrometer mode does not require the full pulse compression. 

For example, a partial compression of the chirped waveform to 20 ns 
(17.3 m surface range resolution) would be quite adequate. The 
excess bandwidth, of course, is still useful for reducing the fading 
variance. With 2i i ns resolution, something like 1024 sample gates 
would adequately sample the return (17.7 km of surface range). The 
spectrometer mode data can be merged with the altimeter mode data 
stream in a way that is compatible with the existing instrument's 
data system. For example, the spectrometer data might consist of 58 
6% bandwidth spectral estimates covering the wavelength range 50-1000 
m output at a nominal 4 frames per second. This data can easily be 
merged with the altimeter data without exceeding the 10 kbs -1 data 
rate of the existing system. Thus, on-board recording, and hence 
fully global coverage is possible. 
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The above figures clearly Indicate that these measurements are feasible. 
Yet, have we missed something, erred. in our thinking somehow? The 
aircraft data to be presented below convince us that we have not. 


4. AIRCRAFT VALIDATION 

4.1 The Fall '78 CV-990 Mission. 

The Goddard K u -band Short-Pulse Radar was built up from the Geos-3 
satellite altimeter breadboard obtained from General Electric Co. in 
1974, and it shares the following characteristics with the space- 
craft instrument: 


Frequency: 
Pulse type: 

Pulse length: 
Peak power: 
PRF: 

Detection: 


13.9 GHz 

Linear FM, 100:1 pulse 
comp ression 

12.5 ns compressed 

2.5 KW 
100 Hz 

Noncoherent square la w 


Prior to 1978, the radar was flown on several aircraft missions with 
fixed-azimuth, variable elevation antennas. A description of the 
Goddard radar as it was configured in 1975 is given by Le Vine et 
al . (1977). A major breakthrough in our program occurred in 1978 when 
we had an opportunity to fly piqqyback, free of charge, on the 1-month 
long, Convair-990 Nimbus-7 Underflight Mission. For this mission, 
one of the fixed-azimuth printed-circuit antennas was modified (by 
sawing it in half) and adapted to an azimuth scan. Also, the data 
system was redesigned to allow continuous recording at the full PRF. 
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Figure 4 shows the rotating antenna installed in the CV-990's 
instrument 'sled'. It is shown surrounded by a cylindrical baffle 
which was designed to protect a neighboring radiometer from possible 
RFI. Also shown in figure 4 is a 12° X 12°, nadir-directed 
rectangular horn antenna, which served In our Instrument's 'altimeter' 
mode. The nadir horn and rotating antennas are shown connected by a 
waveguide switch; this switch could be activated by a mode-change 
conaiand from the radar's control panel in the aircraft cabin. The 
rotating antenna characteristics are: 

Boresight incidence angle, e 0 = 15.3° 

Azimuth beamwidth, 0<j, = 4° 

Elevation beamwidth, 0 q = 10° 

Rotation rate, * = 6 rpm 

The boresight angle was chosen so that an elevation sidelobe at 
15.8° to the main-beam axis would be directed toward nadir. The 
return from the sidelobe, which was recorded in the same frame 
as the main-beam return, allowed us to calculate the range on the 
surface without having to calibrate for absolute time delay. This 
is important in the relatively low-altitude (10 km) aircraft geometry 
where a rather broad elevation beamwidth is required to generate 
a large enough range footprint extent for wave number resolution. 

Thus in the aircraft geometry* wave-front curvature in the elevation 
plane is not negligible, and if not properly accounted for, the 
curvature will result in a considerable dispersion of the surface 
wave number. If t is the time elapsed from the time of the nadir 
sidelobe return, then given the aircraft altitude from the plane's 
operational altimeter, the surface range x as measured from the nadir 
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point can be calculated according to the equation 

x 2 + H 2 * (ct/2 + H) 2 (4.1) 

Of course, in the satellite measurement geometry, we can linearize 
(4.1) to get (2.14). The aircraft measurement geometry is Illustrated 
in figure 5. At the nominal aircraft altitude of 10 km, the footprint 
dimensions are L** * 1500 m and Ly* * 700 m approximately. Because 
of the rapid roll-off of o° with e , the backscattered power peaks 
inward of the boreslght angle. Generally, the peak return occurs in 
the vicinity of 13° incidence. The metallic baffle, and the poor 
radome environment in general spoiled the gain pattern to such an 
extent that we have not attempted to measure o° either as a function 
of elevation or azimuth angle. Figure 6 is an example of the (azi- 
muthally averaged) average backscattered power profile: the large ~ 3° 
ripple near the beam axis caused by diffraction by the baffle is 
obvious. The poor gain pattern is unfortunate as, Ideally, we want 
to estimate the tilt model sensitivity term 8 1 np/a tan© directly 
from the observed cross-section roll off. The gradient of the slope 
pdf and mean square slope are Internal parameters of the measurement; 
yet, in the analysis to follow we shall have to rely on external 
parameters in order to calculate the tilt model sensitivity a . That 
is, we will have to use a mean relationship between the mean square 
slope and the buoy-observed wind speed in order to verify the prediction 
of the tilt sensitivity coefficient (2.11). 

The digital data system consisted of a high-speed waveform sampler 
(Biomation Co.), two 6 k-byte buffers, and a high-speed (75 ips) 
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1600 bpl tape drive. The Blomation sample gates were selectable and 
could be set to 2, 5, 10, or 20 ns. Quantization was 6 bits and the 
maximum frame size was 1024 samples. Generally, we recorded In the 
spectrometer mode at a 5 or 10 ns rate, taking 512 samples at the 
full PRF. Shaft encoder and other housekeeping data were recorded In 
the first two tape tracks. 

The Fall '78 Mission took In 19 flights of approximately 5 hours 
duration In the period October 24 - November 19, 1978. About half of 
these flights were over Ice, the remainder over water. Approximately 
50 2400 ft tapes were written with ocean backscatter data; roughly 
half these data were taken In the Instrument's spectrometer mode, the 
remainder In the altimeter mode. The spectrometer files are h" nd 
large 1-2 minutes long. The 1 minute files, amounting to onl r 
antenna rotations, are a bit short on equivalent DOF and consequently 
the spectra from these files are noisy. 

In this paper we are only concerned with validating the technique, 
and so we shall be examining only a small subs#, of the Fall '78 
Mission data set; that Is, we shall only be examining those files 
for which we have corroborative 'surface truth'. Table 1 is a summary 
of the surface truth data set (spectrometer mode). This data set 
consists of overflights of three types of wave-recording buoys, 
including two N0AA data buoys (N.E. Pacific), a Waverlder (Norwegian 
Sea), and a pitch-roll buoy (N.E. Pacific). Colocation was generally 
within 100 km spatially and within 1 hr temporally. 

4.2 Data analysis 

The digital flight tapes were reformatted and compressed by averaging 


3 consecutive pulses. Also, the spectrometer mode data were 
standardized to 10 ns resolution. Figure 7 is an example of the 
backscatter data contained on the reformatted tapes. The figure 
shows 1500 pulse returns, intensity-coded, and stacked vertically, 
on a CRT display. These (essentially raw) data were further 
processed on a general purpose computer as follows; 

1) The equally spaced array in time W(mAx), m * 1, 2,. ..512, 
is converted to an equally spaced array in surface range according 
to (4.1). The nominal surface range resolution Ax B 8.1 m at 

13° incidence; however, it was more convenient to array the data 
In equally spaced 12 m surface range bins. 

2) The geometrically corrected data are then subjected to two 
algorithms. In the first, no motion correction is applied; the data 
are smoothed in range, and averaged over the several rotations of 
the antenna* This produces an estimate of the average backscattered 
power W 0 (x,4) . In the second alqoritnm, the motion compensation 

is applied. That is, the array is transformed according to 
x + x + Vt cos<t. The data have been processed in 15° azimuth 
blocks; given the 6 rpm rotation rate this means we have integrated 
N « 42 independent pulses. (The pulses are generally independent 
even for forward/aft looks since the sea-Doppler spread is generally 
greater than 100 Hz.) The decision to process in 15° blocks was made 
some time ago; this seemed reasonable as the nominal directional 
resolution is about 17° (200 m wave) and tests indicated no loss of 
signal strength. However, this is about equal to the angle subtended 
by the footprint, and therefore the 15° is at odds with our prior 
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dictum (2.23) to allow no more than half a beamwldth's movement 
in azimuth. Thus, the problem of the appropriate integration time 
for the aircraft data is still somewhat open, and requires a second 
and more careful examination. (As discussed in Section 2, choosing 
too long an integration time will not affect the spectral shape, 
only the signal strength; that is, the signal strength will start to 
drop as the Integration becomes Incoherent with respect to the 
modulation signal . Thus, it is possible that we may have to reexamine 
our conclusions about the precise values of the sensitivity 
coefficient .) 

3) The accumulated N-pulse average Is normalized by the estimate 
of the average power W 0 (x,4>) and unity is subtracted. (The 
azimuthal dependence of W 0 Is probably mostly due to asymmetry 

in the radome environment.) The data are then rewindowed by a cosine- 
squared window. In the high altitude range (8-10 km) the window end 
points are taken to be x « 800 m and x * 3872 m, and at the low 
altitudes (4-6 km) these values are halved. The midpoint of the 
window corresponds roughly to o ■ 13. 5° » 

4) Estimates of Pm(K,$) for each 15° azimuth block are computed 
using a 256 point Fast Fourier Transform. These estimates are then 
ensemble-averaged over the several rotations of the antenna. 

Figures 8a-f are polar contour plots of the processed directional 
spectra P42MK, 4) In units of meters where Pn* s 4tiPn is the one-sided 
(in K) spectrum as a function of wave number in cycles per meter. 

The noise background has not been subtracted in these plots, but as 
the SNR is quite high (+10-20 dB), the background noise level is 
insignificant. Thus, these spectra cat be viewed as directional 
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slope spectra (one need only supply the 'calibration constant' a ). 
Figure 8a Is the directional spectrum of a storm sea* significant 
wave height H s » 9.4 m* dominant wavelength » 330 m. This spectrum 
(Tape 36/FI le 1) was produced from a long run, and hence It Is very 
stable. Figures 8b and e are examples of blmodal spectra. Figure 8c 
represents a fairly low sea state (H s = 1.9 m). The 'rattiness' of 
this spectrum Is characteristic of the noisy spectra from the short 
files (< 6 rotations). Figure 8f Is an Interesting example. It Is 
the observed spectrum of a unidirectional, ~ 330 m monochromatic swell 
running under fairly light winds. Visually, from the vantage point 

of 8.4 km, the sea surface had the striking appearance of a diffraction 

* 

grating. The crest length was for all practical purposes Infinite 
(> 40 km). (The eye can easily detect any deviations from a straight 
line.) It follows then that if the radar system's response were 
perfect, the observed spectrum would be a symmetrical pair of delta 
functions. The blossoming in figure 8f thus represents, In effect, 
the system transfer function (for this 330 m wave).* According to 
(2*2), the directional dispersion should be approximately 30° for 
this wave; indeed, this Is the amount of dispersion seen in the 
spectrum. The wave number dispersion is determined by the finite 
range footprint dimension of 1500 m. If we figure roughly a dispersion 
SK ~ 1/L X * 7 X 10~* cpm, this accords with the observed wave 

number dispersion in figure 8f. 


*Actually, this is the spectral window; one might also refer to 
figure 8f as a map of the spectral 'point spread function', to 
borrow a term from the visible imagery people. 


Lastly, respecting figure 8, we note the asymmetry that is evident 
in several of the spectra. This is obviously related to asymmetry 
in the wave-slope distribution. What is significant here is that the 
asymmetry is by and large rather small, and this would indicate that, 
by and large, the second-order hydrodynamic effects are small. 

Figure 9 is a series of cuts through the directional spectrum Pn(K,<&) 
of Tape 36/File 1. The figure is Intended here to show that the 
foreward face of the (slope) spectrum is sharply defined and stands 
clearly away from the residual antenna pattern energy near dc. This 
rather 'clean' situation will not obtain for the nondi recti onal 
spectra we shall be computing in the following. 

4.3 Absolute nondi recti onal comparisons. 

For these comparisons, we will need to take a closer look at the 

residual fading spectrum which must be subtracted from P^j to give 

the directional modulation spectrum P m . Because of the nonlinear 
time-delay versus surface range relationship (4.1) that obtains in 
the aircraft geometry, the formula (2.18) for P w will not be exact. 

The pulse spectrum in the surface wave number domain will be similar 

to the pulse spectrum only when there is a linear relationship between 
surface range and delay time. Nevertheless, (2.18) may stand as a 
fair approximation. If we assume as a nominal range resolution the 
resolution at 13.5° incidence, i.e., Ax = 8.14 m, then (2 .18 ) gives 
(with N = 42), 


4ttP w = 0.58 [m] exp[-0.5(K/0 .033) 2 ] 


(4 .2) 
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where K Is given in cpm . Figure 10 is a plot of the azimuthal ly- 
integrated value of P^ in the wave number band 0.0218-0.025 cpm 
(center frequency * 0.19 Hz) for several files, where the buoy- 
observed wind speed is used as an ordering parameter. From the plot 
one sees that an extrapolation of the observed Pn to zero wind 
speed (and hence, presumably, to zero modulation depth) yields the 
value of the residual fading variance predicted by (4.2). The pulse 
spectrum roll-off between zero wave number and 0.023 cpm (42 m wave 
length) is about 20%. 

With P ni computed from P^ by subtracting (4,2), the directional 
height-frequency spectrum S(f,4>) is computed using the tilt model 
solution (2.8). Assuming the linear, deep-water dispersion relation 
ship it follows that 


S(f.*) » (2/af)P m (K,*) 

Where S is the polar-symmetric spectrum expressed in m 2 /Hz/rad. In 
computing (4,3), the measured modulation spectrum is symmetrized 
according to 

P m * 0.5 [P„,(M) + P m (M+180»)] 

Symmetrizing the spectrum, of course, has the advantage of doubling 
the DOF. The nondirectional spectrum is now computed according to 

12 

■5(f) - (1/12) l 4»P_(M 1 ) 

1=1 m 1 

= %f m (K)/af 


(4.3) 


(4.4) 


(4.5) 
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As we have pointed out before, the sensitivity coefficient a should 
Ideally be calculated on the basis of the observed cross-section roll- 
off, but the poor antenna gain pattern has made this virtually Impossible, 
Thus, we will have to make do here with an Indirect means of verifying 
the tilt-model predictions. First, let us measure a by taking the 
ratio of the area under the radar spectrum to the area under the 
buoy spectrum. That Is, let 

0.2 Hz 

“mess ' (*•«> 

/?c *(f)df 5 UW 

where f c Is a low-frequency cutoff, and the significant wave height 
H s s 4 <^ 2 >i /2 (where the low-frequency deficit Is understood). The 
low-frequency cutoff Is necessary because the antenna pattern effects 
can be quite severe In the nondlrectlonal spectra computed according 
to (4.5). This Is because the antenna pattern energy Is omnidirectional 
and therefore makes a large contribution to the azimuthal average. 

This should be apparent from figures 8, 9 and 11. (Also, the 1/f 
factor In going to the height spectrum will magnify any errors 
In the specification of the fading noise background level, producing, 

In effect, a 1 1/f noise* background.) 

Table 2 list? the measured alphas for 10 files, representing 
basically 7 Independent observations. The cutoff frequency In each 
case Is chosen arbitrarily as the frequency of the minimum between 
the dc and the spectral peak. Figure 11 gives five examples of the 
Inferred directional height spectra based on the measured alphas. 

Figures 12a«e compare the Inferred nondlrectlonal spectra with buoy 
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observations. The five examples shown are plotted autoscaled, linear- 
linear; they cover a range of sea states from H s » 1 .9 m to 9.4 m, 
and Include a variety of spectral forms. It Is seen that the agree- 
ment Is generally excellent over the entire range from f c to 0.2 Hz, 
The minor discrepancies that are apparent In some of these comp- 
arisons can be attributed, for the most part, to sampling variability, 
geophysical variability (colocation error), and antenna pattern 
contamination . One does not need to look for explanation In terms of 
second-order scattering effects . These effects are by and large so 
small as to be masked by the larger errors . For example, consider 
samDllna varlabllltv. The 90£ confidence interval on the buov soectrum 
In figure 12a is (0.6"?, 1 .9 S) . Thus the confidence Interval on the 
peak of the spectrum is 130% of the full scale of the figure! The 
pattern contamination Is evident in all figures, but It Is only 
severe in the case of figure 8c. In this case, both the frequency and 
the spectral density are low. If one examines all the figures, it 
would appear that the antenna-pattern-related dc component has a 
spectral density of about 2 m 2 /Hz in the vicinity of 0.08 Hz; this 
would account for the apparent discrepancy In figure 8c. Yet, there 
still may be a real low-frequency Whitening due to second-order 
scattering effects (Intermodulation products;— see J); however, 
it Is not possible with the present data to distinguish the real 
whitening from antenna pattern contamination. One would need a 
much larger footprint to reduce the dc component. There does on 
the other hand seem to be evidence of second-order effects on the 
high-frequency side: these are manifest In the slightly more rapid 
roll-off of the radar spectrum. This apparent 'droop 1 Is related 
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to (2.13) and has been predicted by 0. (We will defer for the 
moment a discussion of the colocation errors.) 


From these comparisons, we conclude that the measurements can be 
made with good spectral fidelity. But, how about absolute spectral 
levels? Let us check the measured alphas. How do they compare with 
(2.11)? An Immediate check on the l/Ly dependence is possible with 
files 89/2 and 89/3. The ratio of the altitudes is 9. 3/4. 5 = 2.07; 
the ratio of the alphas is 4.90/2.47 ■ 1.98 (a 4% difference). Thus, 
the asymptotic l/Ly dependence is verified. Now, using (2.11) let us 
infer a mean square slope from the measured alphas. Let us assume a 
nominal incidence angle of 13°. Then, inverting (2.11) we obtain 


<|vc| 2 > = 


2 tan 13°. 



- cot 13° 


(4.7) 


The inferred mean square slope values are tabulated In Table 2 
and plotted in figure 13 as a function of the buoy-observed wind 
speed. No corrections were made for anemometer height or atmos- 
pheric stability. An eyeball regression yields the relationship 

<|Vc| 2 > = 0.0028 U [ms" 1 ] + 0.009 (4.8) 


for the wind speed range of approximately 5 to 20 ms" 1 . Equation 
(4.8) is in perfect agreement with the K u -band scatterometer 
data analyzed by Wentz (1977), at least up to the largest wind 
speed (12 ms" 1 ) in Wentz's data set. In fact, Wentz's data are not 
plotted because they are sensibly no different, and would only 
clutter the plot. Equation (4.8) also agrees with Wilheit's (1979) 
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analysis of passive microwave data. It predicts slope variances 
that are ~ 60% of the optical values reported by Cox and Munk (1954) 
and which, interestingly, lie between their 'clean* and 'slick' 
surface observations. Our inferred mean square slope values (4.8) 
are thus consistent with our knowledge of what the K u -band effective 
slopes should be and strongly support our conclusion that the tilt 
model solution (2.8) is a correct first-order relationship.* 

Now let us use the regression result (4.8) in (2.11) to compute a 
'theoretical' alpha. The 'theoretical' alpha can then be used to 
compute a ' radar-inferred 1 absolute height spectrum and significant 
wave height. Table 2 Is a tabulation and figure 14 is a plot of 
the results for the 'inferred' wave height for the seven independent 
cases analyzed. Over the wave height range 1.9 m to 9.4 m the mean 
difference between the radar-inferred and buoy H s is 0.00 m (sic ) 
and the rms difference is 0.16 m. This is truly remarkable considering 
that a) we are using only a first-order, back-of-the envelope theory, 

b) our measurement geometry is not ideal (broad elevation beamwidth), 

c) we have had to rely on external parameters (buoy wind speeds) 
rather than internal parameters (cross-section roll-off), and lastly, 

d) the data are subject to sampling variability as well as geophysical 
variability (colocation errors). Some information as to the last 
source of error is available to us through the instrument's altimeter 
mode. The altimeter mode algorithm consisted of epoch realignment, 
and an iterative least -squares fitting of an error function to the 
leading edge of the average pulse return. H s was computed from 

the measured temporal dispersion a according to 

*A factor of two error in the computation of P m from the FFT led 
to the wrong result for mean square slope in Jackson ert aK (1981). 
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Hs * [4c2<j2 - ty2]l/2 (4.9) 

where Hp * 4.91 m (compare with Fedor et aU» 1979). The altimeter 
wave heights are shown In the last column of Table 2 and are plotted 
In figure 15. The altimeter H s show a positive correlation with 
the spectrometer-mode minus buoy H s residuals Indicating that 
colocation errors are a significant component of the error budget. 

This probably explains the discrepancy in the blmodal spectrum 
comparison of figure 12e. However, there does remain the possibility 
that the sensitivity coefficient Is not truly isotropic. Data from 
future flights with an Improved radome should enable us to resolve 
this question. 

4.4 Directional comparison. 

The single directional comparison available to us Is with a pitch- 
roll buoy deployed by NOAA from the Canadian vessel Quadra occupying 
Ocean Weather Station PAPA in the Pacific. Five directional spectra 
obtained from the buoy near the time of overflight were provided to 
us by the NOAA Atlantic Oceanographic and Meteorological Laboratories. 
Two radar files, 89/2 and 89/3, are available for comparison, but 89/3 
Is preferred because of Its greater length, even though It Is at the 
lower altitude. The directional height spectrum from 89/3 Is shown In 
figure lid, and the nondi recti onal comparison Is made in figure 12d. 
Because of the basic 180° ambiguity In the radar spectrum, a true 
direction of wave travel must be assumed If first harmonics are to 
be compared. (As suggested by W. J. Pierson (personal communication), 
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probably the most objective way to compare the buoy and radar data 
Is on the basis of the first two angular harmonics since these are 
linearly related to the slope auto- and cross-spectra obtained from 
the buoy record.) To facilitate the comparison* the radar spectrum Is 
set Identically equal to zero In the half space which the buoy 
Indicated to be opposite to the direction of dominant wave travel . 
Bocause the sea Is rather broadly spread In this case (apparently 
due to wind turning) and because there Is swell running at ~ 90° to 
the dominant wave direction, this ad hoc procedure will necessarily 
result In some error, particularly at the low swell frequencies. 
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angular distribution are compared to the buoy data In figure 16. The 
reason that there Is a bias In direction is apparently due to the 
fact that the buoy's compass card had slipped, and true direction 
was figured after the fact from the general meteorological conditions 
(W. Mcleish, personal communication). There Is evidently good 
agreement In first harmonic amplitudes from the peak to 0.2 Hz, but 
the second harmonic amplitudes disagree by almost a factor of two 
over most of the range. Why? Our guess Is that the problem Is with 
the buoy, but the reader may judge for himself. A standard practice 
In the analysis of pitch-roll buoy data Is to assume a model angular 
distribution function such as the following 


D(f,o) = cos 2S [($-$ 0 )/23 


(4.10) 


The first and second harmonics of this distribution can then be 
equated to the observed first and second harmonics to give two 
estimates, si and S 2 , of the spreading parameter s as a function 
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of frequency (Cartwright, 1963). If the two s-parameters are 
Identical* then the actual angular distribution Is perfectly 
represented by the model function (4.10). Figure 17 Is a plot of the 
two s-parameters for the radar data. They are practically equal over 
the entire range of frequencies. This makes sense since,, according 
to figure 10* the radar data do appear to be consistent with a cosine- 
power distribution. But now, the buoy data do not appear to be 
consistent with a cosine-power model . This Is apparent from figure 
19, which shows the radar data falling neatly on Cartwright's (1963) 
curve while the buoy data scatter over a large area above the curve. 

This consistency of the radar data versus that of the buoy data 
naturally Inclines us to believe more In the radar results, 
especially as we are biased In that direction to begin with. 
Nonetheless, we should avoid drawing hasty conclusions. In truth, 
the disagreement here is not as great as figure 19 might lead one to 
believe. First, recall that there is agreement in the first harmonic 
and si values. At the peak, si *» 4.9, which Implies a half-power 
spread of 85° . Second, the half-power spread Is not very sensitive 
to the s-value, at least for s > 5 or so. While the buoy s-values 
scatter considerably, they average about sg = 10 near the peak. 

This Implies a spread of 59° which, really. Is not that far away 
from 85°. Further, the average of the two values, 72°, is quite 
close to the observed radar half-power spread of 67° at 0.114 Hz. 

Still, we are skeptical of obtaining from buoys the kind of 
Intercomparison data we need to gain a full understanding of the 
limits of accuracy of these radar measurements, especially as 
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concerns the radar's ability to measure directional spectra of 
complex seaways. Another new aircraft radar development, the NASA 
Surface Contour Radar (Walsh et al_., 1981), promises to be able to 
provide exactly the kind of high-resolution Intercomparison data we 
now require at this point. 

5. CONCLUSION 

We have described a rather simple microwave radar technique for 
measuring directional wave spectra. We have shown that satellite 
measurements on a truly global scale are possible with this technique 
nlque, and we have, In our opinion, provided a firm theoretical 
and experimental basis for the technique. While further aircraft 
experimentation 1$ warranted and desirable we are going to be 
limited by the constraints of the relatively low aircraft altitudes. 
That Is, at aircraft altitudes in a near-nadir geometry, there is 
no way to avoid both a relatively broad beamwldth on the one hand 
and a relatively small footprint on the other. Both of these factors 
make it difficult to investigate second-order scattering effects: we 
may have the theory, but we will be unable to satisfactorily verify 
it. Without a space experiment, there Is little reason to develop 
the theory of measurement any further; that Is, to go beyond Jackson's 
(1981) theory. 

We do Intend however to improve the aircraft radar system. For 
example, we are replacing the 12° nadir horn antenna with a much 
broader-beamwldth standard gain horn. This should permit accurate 
measurements of the Isotropic cross-section roll-off and mean 
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square slope for better estimates of the first-order sensitivity 
coefficient. This type of measurement has been described by 
Hammond et_ a]_. (1977). 
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TABLES 

Table 1. Fall *78 Mission Surface Truth Data Summary (Spectrometer Mode). 
Table 2. Measured Versus 'Theoretical* Sensitivity Coefficient. 
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Figure 1 
Figure 2 

Figure 3 
Figure 4 


Figure 5 
Figure 6 

Figure 7 


Figure 8 


Figure 9 


FIGURE CAPTIONS 


. a) Typical satellite measurement geometry, b) Scan pattern on 
the ocean surface, 700 km altitude, 10° Incidence angle, and 
3 rpm scan rate. 

. Illustrating directional selectivity by phase-front matching 
of em and ocean wave components. For the rectangular Illumi- 
nation pattern Illustrated here the angle of the first null 
Is as Indicated. 

. Simple tilt model of reflectivity modulation, 

. Antennas mounted In the CV-990 instrument "sled". View Is 
upward, looking Into the sled with the radome cover removed. 

The rotating antenna Is surrounded by a cylindrical baffle. 

The nadir-pointing horn antenna is connected to the rotary 
antenna by a wave guide switch. The other antennas shown 
belong to the SMMR Simulator microwave radiometer. 

, Aircraft measurement geometry a) elevation view, b) plan view. 

. Azimuthal ly-averaged, average backscattered power profile, 

Tape 37/File 1. The upper and lower dashed curves represent 
respectively the average maximum and minimum values over 360° 
of azimuth. 

. CRT display of 3-pulse average backscatter data. The sample 
gate setting Is 10 ns and the frame size is 512. The display 
represents 512 consecutive 3-pulse averages stacked vertically. 
The tic marks are placed every 250 pulses, or, equivalently, 
every 90° of antenna rotation. The "S" pattern is the result 
of antenna rotation combined with aircraft motion. 

. Polar contour plots of the (one-sided) directional spectra 
4 ttP 42 (K,*) In units of meters. Top of the figures Is zero 
degrees radar azimuth « aircraft heading. The wave number 
rings are spaced every 0.005 cpm (first ring, 200 m wavelength, 
second 100 m, and third 66.7 m). The labeled contour level 
(CLEV) Increments are low by a factor of 2; they should read 
for a) Tape/File « 36/1, CLEV = 9.08 m, b) Tape/File « 91/6, 

CLEV « 5.00 m, c) Tape/File « 85/10, CLEV = 3.00 m, d) Tape/File 
* 89/3, CLEV * 5.16 m, e) Tape/File * 94/2, CLEV * 3.76 m, 
and f) Tape/File = 90/7, CLEV = 2.8 m. 

. Radial cuts through the directional spectrum of Figure 8a. 
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Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 

Figure 


10* Estimation of the background fading noise level . The level 
of the azlmuthally averaged directional spectrum In the wave 
number band 0.0218-0.0250 cpm for 12 files Is plotted versus 
buoy-observed wind speed. The data from the adjacent files 
45/2 & 3, 86/4 4 6, 89/2 4 3, and 90/7 4 9 have been averaged 
together. An eyeball extrapolation to zero wind speed yields 
the value of the residual fading variance spectrum predicted 
by equation 4.2 

11. Polar -symmetric directional height spectra S(f,$) corresponding 
to figures 8a-e, The frequency rings are 0.1 H z and 0.2 Hz. 

The equally-spaced contour levels (CLEV), In units m 2 /Hz/rad, 
have been multiplied by a factor of 2ir. 

12. Comparisons of radar-inferred (solid line) and buoy (circles) 
nondl recti onal height spectra S(f) using the measured alphas, 
a-e) correspond to a-e) of figure 11. In d) the circles and 
squares stand for the buoy records at 2138 Z and 2104 Z respect- 
ively. The 2138Z data have been smoothed by a 2-polnt average 

In the vicinity of the peak. 

13. Mean-square slopes Inferred from the measured alphas versus 
wind speed, 

14. Radar-Inferred versus buoy significant wave height (spectrometer 
mode). Data from adjacent files have been averaged together; 
cf. Table 2. 

15. Radar-Inferred versus buoy significant wave height (altimeter 
mode ) . 

16. First and second harmonics of radar angular distribution 

compared to pitch-roll buoy data. Solid line ( _)1s from 

radar file 89/3; dashed ( ) is from average of five buoy 

records; and dotted (••••) Is from buoy record at 2104 Z. 
a,b) phases; c,d) amplitudes. 

17. Radar-observed spreading parameters corresponding to the 
first and second harmonics shown In figure 16. 

18. Radar-observed directional distributions (Tape/File * 89/3) 
compared to cosine-power distribution. The symmetrized 
radar distribution Is set identically equal to zero between 
0° and 180°. 

19. Radar (bullets) and buoy (cpen circles) first and second 
angular harmonic amplitudes compared with theoretical relation- 
ship of Cartwright (1963). Hie data, from radar file 89/3 

and buoy record 2104Z, are contained In the frequency range 
0 .078 Hz to 0 .145 Hz . 
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